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1. Introduction 

The shell map is a very simple representation of the structure of foams, com- 
bining the geometrical (random tiling) and dynamical (loss of information 
from an arbitrary cell out) aspects of disorder. We will illustrate it and give 
several examples, including a few arising from discussions in Cargese. This 
chapter is written by following the main lines of two previously published 
papers [1, 2]. 

In Nature, space-filling disordered patterns and cellular structures are 
widespread [3, 4]. These structures (froths) are partitions of D-dimensional 
space by convex cells. Disorder imposes that each vertex has minimal num- 
ber of incident edges, faces and cells (D + 1 edges incident on a vertex, D 
faces incident on a edge, D — 1 cells incident on a face, in D-dimensions, 
Fig.l). In this respect, a froth is a regular graph, but the number of edges 
bounding each face, the number of faces bounding a polyhedral cell, etc., are 
random variables [5]. Minimal incidences implies also that the topological 
dual of a froth is a triangulation (Fig.2), a useful representation of packings. 
Indeed, for any given packing, or point set, one can construct the Voronoi 
tessellation [6], which is a space-filling assembly of polyhedral cells. When 
the starting points are disordered (no special symmetries) the Voronoi tes- 
sellation is a froth. The space filled by the froth can be curved (Fig.l). This 
is for instance the case in amphiphilic membranes [7], fullerenes, the basal 
layer of the epidermis of mammals [8] or the ideal structure of amorphous 
materials [9, 10]. Disorder does not necessarily imply inhomogeneity. On the 
contrary, in many cases, disordered froths are very homogeneous (cells with 
very similar sizes and regular shapes). This is -for instance- the case in the 
epidermis, where the biological cells have homogeneous sizes and isotropic 
shapes but the structure is disordered. Indeed, the disordered arrangement 
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Figure 1. An example of two dimensional disordered cellular structure (froth) generated 
by cell division and coalescence transformations (left). A froth on a curved space (right). 

is the one which best guarantee both the partition of the curved space into 
similar cells and the invariance under mitosis and detachment. 

The interplay between disorder and curvature is illustrated in this chap- 
ter by representing the froth as organized in concentric layers of cells around 
an arbitrary central cell [1] (Figs. 3 and 5). The structure is built from the 
central cell outward like an ever expanding jigsaw puzzle without boundary. 
The radial map, from one spherical layer of cells to the next, is the logistic 
map [11], and the geometrical tiling is expressed mathematically as a dy- 
namical system [12]. The isotropy of the disordered structure is expressed 
locally by averaging over each layer. The over-all translational invariance is 
manifest in the independence of the structure and properties on the choice 
of the central cell. The radial map from one layer to the next includes both 
effects of disorder and of space curvature. 

In a two dimensional froth, given a cell with n edges, one can define its 
topological charge as q = 6—n. The total charge of a froth with N cells is the 
sum over the charges of each cell: Qt = J2i Qi = (6— {n))N = 6%, with \ the 
Euler-Poincare characteristic of the manifold tiled by the froth [13, 14]. The 
charge is a topological invariant, it cannot be generated or destroyed and 
the local topological transformations in the froth redistribute it between 
adjacent cells. The total charge is equal to zero in Euclidean froths, it is 12 
for froths on the surface of a sphere and it is negative in spaces with negative 
Gaussian curvature. In two dimensions, it is therefore possible to define the 
curvature of a surface tiled by a froth by analyzing the local topological 
configuration of its tiles. In three dimensions, for N — > oo or in closed 
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Figure 2. Froths and triangulations are dual structures. Triangulations are useful rep- 
resentations of packings. 

elliptic froths, the Euler's relation is homogeneous (% = 0, independently 
on the space-curvature) and from the Gauss-Bonnet formula [14] it is not 
possible to distinguish the global curvature of a tiled manifold from the 
local average properties of its tiles. On the other hand, we show that, also 
in three dimensions, the map makes possible to define the curvature of 
the space from local topological configurations [1]. This is done by simply 
computing the number of cells in successive layers. 

In this chapter, we describe the map which gives the number of cells 
Kj in a layer distant j from the central cell as function of the average 
topological properties of the cells in the previous layers (section 2) [1]. We 
discuss the link between the map and the space curvature giving examples 
in two and three dimensions (section 3) [1, 2]. By using the map we exploit 
the freedom of constructing a froth with different local topological config- 
urations and determine the three dimensional Euclidean structures which 
maximize such freedom (section 4) [1, 2]. 

2. From a cell to the whole froth, a topological map 

All froths can be studied as structured in concentric layers of cells which are 
at the same topological distance from a given central cell. The topological 
distance between two cells is the minimum number of edges crossed by a 
path from one cell to the other [1]. The layers are closed rings of irregular 
polygons in two dimensions and spherical caps of irregular polyhedra in 
three dimensions. The cells making the layer j can be distinguished in two 
categories. Some cells have simultaneously neighbours in the layers j — 1 
and j + 1, these cells make themselves closed layers and constitute the 
skeleton of the shell-structure. Other cells (or clusters of cells) are local 
inclusions (topological defects) between the layers of the shell-skeleton (they 
don't have neighbours in the layer j + l). The shell-skeleton is itself a space- 
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Figure 3. An example of SSI froth. The number of cells in each layer can be calculated 
in term of a simple map 

filling froth, hierarchically organized around the germ cell. Once the germ 
cell is chosen, the shell-structure and its skeleton are univocally defined, but 
different germ cells generate different skeletons. We call shell- structured- 
inflatable (SSI) a froth free of topological defects. In this case, the shell- 
structure and its skeleton coincide. An example of a two dimensional SSI 
disordered cellular structure is given in Fig. (3). A three-dimensional regular 
SSI structure (the Kelvin froth [15]) is given in Fig. (4). In this paragraph, 
we first study SSI froths and then generalize the results to generic froths. 

2.1. A TOPOLOGICAL MAP FOR SSI FROTHS 

For SSI structures, a recursive equation gives the number of cells Kj in a 
layer at a distance j from the germ cell, as a function of the number of cells 
in the previous two layers. In two dimensions, one can easely verify that 
[1, 12, 16, 17] (see Fig.3) 

Kj+^SjKj-Kj.! (forj>l), (1) 

where the inflation parameter Sj is related to the average topological prop- 
erties of the cells in layer j, 

Sj = (n)j - 4 (2) 

with (n)j the average number of neighbours per cell in layer j. The map 
starts with Kq = and K\ = n = number of neighbours of the central cell. 
For example in the hexagonal lattice one has n = 6, (n)j = 6 and therefore 
sj = 2. In this case the map (1) correctly gives Kj = 6j. 
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Figure 4- (a) An example of 3D cellular system: the "Kelvin froth". A part of its 
shell-network is shown in (b), it is the 2D elliptic set of facets which are bounding 
externally a layer, (c) Construction of the 3D "Kelvin cell" from the 2D shell network, 
(d) The "twisted Kelvin cell" . 



In three dimensions, one can write a recursive equation similar to Eq.(l) 
[1, 2] 

: i • '•, (forj>l). (3) 

with 



sj = i((f)j-G)((n N h-4)-bj-l 



(n N ) 



b s = 



/i-i 

2( 6 .- s . + 6 ((n Ar ) J -4) + l 



(4) 



Here (f)j is the average number of neighbours of the cells in layer j. The 
quantity (n N )j is the average number of edges per face in the 2D elliptic 
set of facets which are bounding externally layer j (the shell-network, see 
Fig.4(a,b) and [1, 2] for details). The map starts with Kq = and K\ = f = 
number of neighbours of the central cell. For example in the Kelvin structure 
one has / = 14, = 14 and {n N )j = 5 + (24j 2 + 24 j + 7) _1 . By using 
the map (3) we obtain K2 = 50, K3 = 110 and in general Kj = 12j 2 + 2, 
which is the correct answer [18]. 

Equations (1) and (3) take a much simpler form if the average topolog- 
ical properties of the cells in layer j are independent on the layer number 
(i.e. (n)j = (n), (f)j = (/) and (n N )j = (n N )). In general, these quantities 
can vary as one goes from one layer to the next (as for the Kelvin case). 
However, in the limit j — > 00 they must converge towards the averages over 
the whole structure (in the Kelvin case (n N )j — > 5). Moreover, since the 
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choice of the central cell is arbitrary, in disordered systems this asymptotic 
behaviour is reached much faster if one averages over the central cell first. 
Finally, the relations (1) and (3) are also valid when the "central cell" is a 
cluster or even an infinite set of cells. For instance the starting configura- 
tion (the cells with j = 0) can be an infinite linear strip of adjacent faces in 
2D or a planar layer of adjacent bubbles in 3D. In this case the asymptotic 
behaviour and the average over the system can be set from the beginning. 

When the average topological properties of the cells in layer j are inde- 
pendent on the layer number, Eqs.(l) and (3), give [1, 12] 




This equation is the logistic map [11]. In two dimensions, Kj = Kj is the 
number of polygonal cells which are making the layer j. In three dimensions, 
for s^2, Kj is the number of polyhedral cells making the layer j plus the 
additive constant c/(2 — s). When s = 2 Eq.(5) cannot be used and Eq.(3) 
must be used instead. 

2.2. TOPOLOGICAL DEFECTS AND NON-SSI FROTHS 

In almost all disordered froths, non-SSI inclusions (topological defects) are 
present. An example of non-SSI froth is given in Fig.5. A defect in layer 
j is a cell (or a cluster of cells) which have no neighbours in layer j + 1. 
Examples of 2D topological defects are given in Fig.5. The fraction (5) of 
defects in layer j respect to the total number of cells, tends to a constant 
as j tends to infinity. Typically, in 2D physical froths one finds S c± 0.1 [19], 
whereas very disordered computer simulated froths have 0.15 < 5 < 0.6 
[20]. 

Maps (1) and (3) apply to SSI froths only. But, for any froth, given 
a central cell, one can always single out the topological defects. Then the 
reduction to the shell-skeleton is made by eliminating the defects. The 
shell-skeleton is itself a space-filling froth hierarchically organized around 
the germ cell, and it is SSI by construction. Maps (1) and (3) are therefore 
applicable on this structure. 

3. Space curvature from the map 

Let start from a given cell of the froth and count its number of neighbours 
Kj at a distance j. We expect different behaviours for Kj v.s. the distance j 
depending on the curvature of the manifold tiled by the froth, as schemat- 
ically shown in Fig. 6. Mathematically, these different behaviours can be 
easely studied from Eq.(5) when Sj = s independent of j. In this case, one 
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Figure 5. Example of non SSI froths (left). Defects, in layer j, are cells which have 
no neighbours in layer j + 1. They are topological inclusions between the layers that can 
be considered as vertex- or edge-renormalizations (right). The topological map doesn't 
applies to such configurations, but it applies on the shell-skeleton, which is the SSI 
structure obtained by eliminating the defects. 



has the following solution of the map (5) [1], 

Kj oc (exp(ipj) - exp(-(pj)j . (6) 

The parameter s separates the solutions into classes. The region |s| < 2 
is associated with tp = icos -1 (s/2). It has bounded, finite trajectories 
in the (J,Kj) parameter space. This is the characteristic behaviour of 
cellular tilings of elliptic manifolds (see Fig. (6c)). For \s\ > 2 one has 
tp = cosh (s/2), and unbounded exponential trajectories (typical of hy- 
perbolic tilings, Fig. (6a)). The point s = 2 separates the elliptic from the 
hyperbolic regions. It corresponds to Euclidean tilings and the solution of 
Eqs.(l) and (3) is Kj oc j D ~ l , with D the space-dimension (this is the 
behaviour expected from elementary geometry, Fig. (6b)). Note that this 
connection between the space curvature and the inflation parameter s ap- 
plies both in two and three dimensions. 

In two dimensions, equation (6) is the topological version of the relation 
(attributed to Gauss) between the circumference C of a circle and its ra- 
dius p measured on a manifold of constant curvature 7Z (where 7Z = Ricci 
curvature scalar = 2 x Gaussian curvature), 



C = -£-(«p(J-f)-«p(-J-|)) . (7) 




Figure 6. The curvature of the space tiled by the froth is obtained from the behaviour 
of the number of cells (Kj) in successive layers as function of the topological distance 
(j) from the central cell. The flat, Euclidean 2D space corresponds to a linear increment 
of Kj with j (b). In elliptic spaces, Kj first increases, reaches a maximum and then 
decreases (c). In hyperbolic tilings this number increases exponentially (a). 

In our case, C represents the topological length (Kj) of the boundary 

of the cluster with topological radius j and is the parameter <p of 

Eq.(6), given here in terms of the curvature. 

Here are some examples of regular 2D froths with different s. To s = 
— I corresponds an elliptic froth made with four triangles, i.e. the surface 
of a tetrahedron. The value s = corresponds to an elliptic froth made 
with six squares, the surface of a cube. To s = 1 is associated an elliptic 
froth with pentagonal faces, the surface of a dodecahedron. It is known 
that the hexagonal lattice (s = 2) is the only regular froth which tiles 
the Euclidean plane. Froths with heptagonal (s = 3) or octagonal faces 
(s = 4) tile hyperbolic surfaces. Two-dimensional disordered tilings can be 
generated by Vorono'i [6] tessellations around points at random on a surface 
with arbitrary curvature. Here the map must be studied experimentally case 
by case. In the Euclidean plane, the linear behaviour of Kj v.s. j predicted 
by our map, has been experimentally confirmed by many measurements 
on natural and computer generated froths [20, 19, 21]. The behaviour of 
Kj in random Vorono'i networks constructed on curved surfaces is under 
investigation [22]. 
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In two dimensions, the classification by the map (5) is identical to that 
provided by the combination of the Gauss-Bonnet theorem [14] and Euler's 
equation: 

J|da = 2vr X = |Q T = |(6-(n))iV = |(2- S )iV , (8) 

where N is the total number of cells in the system. Equation (8) shows 
clearly that the parameter s < 2 corresponds to global positive Gaussian 
curvature, s > 2 to negative curvature and s = 2 to the Euclidean limit. 
In 2D, the determination of the curvature by the map (5) or by the Gauss- 
Bonnet theorem are identical, but the two methods are independent. In 
three dimensions, the Euler's relation is homogeneous and from the Gauss- 
Bonnet formula it is not possible to distinguish the global curvature of 
a tiled manifold from the local average properties of its tiles. However, 
the map (5) is still applicable. It provides therefore a general means of 
describing the curvature from topological considerations, also in 3D. 

Here are two example of 3D froths with different inflation parameter 
s. The regular froth made by packing dodecahedras has (/) = 12 and 
(n N ) = 5, which substituted into Eq.(4) leads to s = 1 , indicating therefore 
that this froth tiles a positively curved space. Indeed, it is known that it is 
a closed structure made with 120 dodecahedras. The Kelvin froth of Fig.4 
has (/) = 14 and asymptotically (n N ) = 5, corresponding to s = 2, which 
correctly indicate that the Kelvin's cells fill the Euclidean 3D space. 

4. Construction of three dimensional disordered structures from 
the map 

Three-dimensional space-filling cellular systems are highly correlated struc- 
tures. Each cell has a different shape and there exist in general very few 
configurations where all cells pack together and fill space without gaps or 
overlaps. The space-filling condition is a very important constraint which 
strongly determines the properties of disordered structures. 

In froths, the incidence numbers (number of edges on a vertex, of faces 
on a edge, etc.) are fixed at their minimal value. By contrast, the coordina- 
tion number (number of faces of a 3D cell, number of edges of a 2D face) is 
a random variable and its average characterize the structure topologically 
[5]. In two dimensions, the Euler relation fixes the average number of edges 
per face (see Eq.(8)). In three dimensions, the Euler relation associates the 
average number of edges per face (n) with the average number of faces per 
cell (/> 
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Figure 7. Three dimensional froths generated by using the map starting from the shell 
network. The shell network is always the result of a superposition of two three connected 
networks, in the examples (a), (d) it is the superposition of two hexagonal lattices. 

So that the two average coordination numbers are related, and one of these 
two quantities is not constrained. 

The representation of the cellular system in concentric layer around a 
given central germ cell can be also regarded as an iterative way of construct- 
ing disordered packings. In the Euclidean space (and in the asymptotic limit 
when (f)j and (n N )j are independent of j), the inflation parameter of the 
map is fixed at s = 2. For this value of s Eq.(4) gives [1] 

<» = 6+ (^Fi ' (10) 

This equation relates the average coordination ((n N )) of a two dimensional 
structure (the shell-network) with the average coordination ((/)) of the 
three dimensional cellular system. The shell-network is the result of a su- 
perposition of two elliptic 3-connected networks, the "incoming" and the 
"outgoing" froths [1]. The pattern of edges constituting the shell- network 
sets the value of (n N ). Therefore, Eq.(10) allows us to construct systemat- 
ically 3D Euclidean SSI froths starting from 2D shell- networks. 

The simplest 2D froth is the hexagonal lattice. The examples displayed 
in Fig.4 and Fig.7(a-e) illustrate the construction of ordered, monotiled 3D 
froths from a shell-network generated by different superpositions of two 
hexagonal lattices (see also [1, 23]). The cell in Fig.4(c) is topologically 
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Figure 8. The shell network for the Z and A15 phases (here the shells are planes). 
In the Z and A15 shell networks there are respectively 12 and 14 pentagons per each 
hexagon. Therefore, (n N ) result 66/12 and 76/15 for Z and A15 respectively. In the Z 
phase, successive shells have alternating configurations (j) and (j + 1). 

equivalent to Kelvin's a-tetrakaidecahedron (it builds up the Kelvin froth 
shown in Fig. 4(a)), and the cell 4 (d), to its twisted variant [24, 25]. Both 
structures have (/) = 14 and asymptotically (n N ) = 5, and fill the Eu- 
clidean space. In Fig. 7(a) is shown part of a shell- networks with 5-sided 
faces, generated by the superposition of two "squeezed" hexagonal lattices 
and its unit cell (b) (see also [1, 23]). This structure has again (n N ) = 5, 
(/) = 14 and is an Euclidean (s = 2) space-filler. With another inter- 
section of the two "squeezed" hexagonal lattices, one generates the shell- 
network shown in Fig.7 (c), it has (n N ) = 4.8 a 3D unit cell (7 (d)) with 
(/) = 16 (8 quadrilaterals, 6 hexagons and 2 octagons). This unit cell is 
a monotile Euclidean space-filler that, as far as we know, was reported 
for the first time in [1]. Recently, this tiling has been proposed [26] as 
the structure of the ternary crystal ThCr^Si^. Fig.7 (e),(f) show an ex- 
ample of an Euclidean shell-structured-inflatable froth made of two differ- 
ent cells. The shell- network in Fig.7 (e) has also two different tiles with 
(n N ) = 32/6 = 5.33... . The associated 3D unit cell (7 (f)) has (/) = 12 
and (n) = 5. This structure is an Euclidean space-filler with s = 2. One can 
in general show that any Euclidean shell-structured-inflatable froth made 
with topologically identical cells can be constructed from a shell-network 
generated by the superposition of two hexagonal lattices. 

We can restrict our attention to a special class of structures where the 
packed polyhedras have only pentagonal and hexagonal faces. In this case, 
the range of variability of (/) is restricted between 12 and 14 by Eq.(9) 
and (10). The lower value corresponds to a structure with only pentagonal 
faces: a packing of dodecahedras, and it cannot be realized in Euclidean 
spaces (it is an elliptic system of 120 dodecahedras). Whether the upper 
value is realizable or not in Euclidean spaces is still an open question. 
(The only known structure with (/) = 14 made only with hexagon and 
pentagons is the Goldberg froth [27], but this packing is hyperbolic [1]). Two 
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Figure 9. Full lines: average number of edges per face in the whole froth ((n)) and in 
the shell-network {(n N )) as function of the coordination number ((/)). The circles and 
diamonds are some examples of known three dimensional cellular structures made by 
packing polyhedra with hexagonal and pentagonal faces only. Several of them (diamonds) 
occurs in nature (the TCP phases of metallic compounds). 



structures with coordinations number equal to 13.2 and 13.333... have been 
obtained by decurving iteratively the curved structure of 120 dodecahedras 
by introducing hexagonal faces [9, 28]. A natural class of ordered structures 
made only with pentagons and hexagons are the 24 known TCP. They are 
Euclidean and have 13.333... < (/) < 13.5 [30, 31, 32, 33, 34]. The TCP 
A15, with (/) = 13.5, is the known packing of polyhedral cells of equal 
volumes with minimal interfacial area per cell [35]. In a previous paper [1], 
we proved that all TCP structures could be described topologically by the 
map (5). In figure 8 are given two examples of shell networks generated by 
combinations of pentagonal and hexagonal tiles only. From these networks 
the structure of the Z phase and of the A15 phase can be constructed. They 
have coordination (n N ) = 66/13 = 5.076... (Z) and (n N ) = 76/15 = 5.066... 
(A15), that substituted into Eq.(10) give correctly (/) = 13.42... for Z and 
(/) = 13.5 for A15. 

Figure 9 summarizes in a plot some of the previous considerations and 
results. The two lines in Fig. (9) are plots of Eq.(9) and (10). As one can 
clearly see, in a given structure, (n) and (n N ) are different in general. The 
packing is organized so that the faces of the shell-network have an average 



13 



number of sides ((n N )) different from that of the whole structure ((n)). But 
there is a point [1] where these two values coincide: 

(n N ) = (n) = (n)* = (11) 

and, from (9) or (10) 

(/)* = 8 + 2V7= 13.29... . (12) 

When (/) = (/)* = 13.29..., an arbitrary cell has the freedom to 
adhere to a layer by any subset of its faces, without adjustment. This 
freedom gives therefore many more local possibilities for building up the 
structure. A coordination number (/) = 13.29... characterizes therefore 
an hypothetical SSI topological packing which fills the Euclidean three- 
dimensional space with the largest possible number of local arrangements, 
that is with maximum entropy. Coordination 13.29... lies in between the 
two froth obtained by decurving dodecahedras and it is inside the interval 
where (/) is ranging in natural structures and computer simulated froths 
[25, 29, 30, 31, 32, 34, 36, 37]. 

Note that, this structure has a coordination number which is irrational. 
It cannot be therefore a periodic crystalline structure. It can only be ap- 
proximated by disordered froths or quasi-crystalline structures. 

The author acknowledge many discussions with N. Rivier. A special thanks to M. O'Keeffe, 
who sent us the preprints of his latests works. This work was partially supported by EU, 
HCM Program, "FOAMPHYS" network, contract ERBCHRXCT940542 and by TMR 
contract ERBFMBICT950380. 
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